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Abstract 



A direct search has been made for magnetic monopoles produced in e + p collisions at 
a centre of mass energy of 300 GeV at HERA. The beam pipe surrounding the interaction 
region in 1995-1997 was investigated using a SQUID magnetometer to look for stopped 
magnetic monopoles. During this time an integrated luminosity of 62 pb _1 was delivered. 
No magnetic monopoles were observed and charge and mass dependent upper limits on the 
e + p production cross section are set. 
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1 Introduction 



One of the outstanding issues in modern physics is the question of the existence of magnetic 
monopoles. Dirac showed that their existence leads naturally to an explanation of electric charge 
quantisation [1]. Magnetic monopoles are also predicted from field theories which unify the 
fundamental forces [2-5]. Furthermore, the formation of a monopole condensate provides a 
possible mechanism for quark confinement [6]. Nevertheless, despite a large number of searches 
[7, 8] using a variety of experimental techniques, no reproducible evidence has been found to 
support the existence of monopoles. Searches for magnetic monopoles produced in high energy 
particle collisions have been made in pp [9-11] and e + e~ [12-17] interactions. This paper 
describes the first search for monopoles produced in high energy e + p collisions. 

The quantisation of the angular momentum of a system of an electron with electric charge e 
and a monopole with magnetic charge g leads to Dirac's celebrated charge quantisation condi- 
tion eg = nHc/2, where H is Planck's constant divided by 2n, c is the speed of light and n is an 
integer [1]. Within this approach, taking n — 1 sets the theoretical minimum magnetic charge 
which can be possessed by a particle (known as the Dirac magnetic charge, g D ). However, if the 
elementary electric charge is considered to be held by the down quark then the minimum value 
of this fundamental magnetic charge will be three times larger. The value of the fundamental 
magnetic charge could be even higher since the application of the Dirac argument to a particle 
possessing both electric and magnetic charge (a so-called dyon [18, 19]) restricts the values of 
n to be even [18]. 

Monopoles are also features of current unification theories such as string theory [2, 3] and 
Supersymmetric Grand Unified Theories [4,5]. Both of these approaches tend to predict heavy 
primordial monopoles with mass values in excess of 10 15 GeV. However, in some Grand Unified 
scenarios values of monopole mass as low as 10 4 GeV [20-22] are allowed. Light monopoles 
are also predicted in other approaches [23-26] and postulates on values of the classical radius 
of a monopole lead to estimates of mass of (9(10) GeV [8]. 

Since the value of the coupling constant of a photon to a monopole (a rn 34n 2 ) is substan- 
tially larger than for a photon-electron interaction (a e w 1/137) perturbative field theory cannot 
be reliably used to calculate the rates of processes involving monopoles. The large coupling also 
implies that ionisation energy losses will be typically several orders of magnitude greater for 
monopoles than for minimum ionising electrically charged particles [27-29]. 

Direct experimental searches using a variety of tracking devices to detect the passage of 
highly ionising particles with monopole properties have been made [7]. Direct searches have 
been made for monopoles in cosmic rays [30] and for monopoles which stop in matter such 
as at accelerators [10] and in lunar rock [31-33]. One method of detection is the search for 
the induction of a persistent current within a superconducting loop [31], the approach adopted 
here. Measurements of multi-photon production in collider experiments [11, 17] allow indirect 
searches for monopoles to be made, the interpretation of which is believed to be difficult [34,35] . 

2 The Experimental Method 

For the direct search reported here we use the fact that heavily ionising magnetic monopoles 
produced in e + p collisions may stop in the beam pipe surrounding the HI interaction point at 
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HERA. The binding energy of monopoles which stop in the material of the pipe (aluminium 
in the years 1995-1997) is expected to be large [36] and so they should remain permanently 
trapped provided that they are stable. The beam pipe was cut into long thin strips which were 
each passed through a superconducting coil coupled to a Superconducting Quantum Mechanical 
Interference Device (SQUID). Fig. [l] shows a schematic diagram illustrating the principle of 
the method used. Trapped magnetic monopoles in a strip will cause a persistent current to 
be induced in the superconducting coil by the magnetic field of the monopole, after complete 
passage of the strip through the coil. In contrast, the induced currents from the magnetic fields 
of the ubiquitous permanent magnetic dipole moments in the material, which can be pictured as 
a series of equal and opposite magnetic charges, cancel so that the current due to dipoles returns 
to zero after passage of the strip. 

The aluminium beam pipe used in 1995-1997 was exposed to a luminosity of 62 ± lpb -1 . 
The beam pipe around the interaction point had a diameter of 9.0 cm and thickness 1 .7 mm in 
the range —0.3 < z < 0.5 m 1 and a diameter of 11.0 cm and thickness 2 mm in the range 
0.5 < z < 2.0 m. During HERA operations it was immersed in a 1.15 T solenoidal magnetic 
field which was directed parallel to the beam pipe, along the +z direction. This length of the 
pipe, covering —0.3 < z < +2.0 m, was cut into 45 longitudinal strips each of length on 
average of 573 mm (~ 2 mm was lost at each cut). The central region (—0.3 < z < 0.3 m) 
was cut into 15 long strips of width ~ 18 mm, two of which were further divided into 32 short 
segments varying in length from 1 to 10 cm. The downstream region (0.3 < z < 2.0 m) was 
divided into 3 longitudinal sections each of which was cut into 10 long strips of width ~ 32 mm. 
The long strips and short segments were each passed along the axis of the 2G Enterprises type 
760 magnetometer [37] at the Southampton Oceanography Centre, UK. This is a warm bore 
device with high sensitivity and low noise level. The samples of strips and segments were 
passed through the magnetometer in steps, pausing after each step, after which the current in 
the superconducting loop was measured. The residual persistent current after the complete 
traversal of a sample through the loop was measured by taking the difference in the measured 
current after and before passage. The readings for each sample were repeated several times. 
This allowed the reproducibility of the results to be studied so that random flux jumps and base 
line drifts could be identified. Any real monopole trapped in the pipe would give a consistent 
and reproducible current step. 

A long, thin solenoid, wound with copper wire on a cylindrical copper former, was used to 
assess the sensitivity of the SQUID magnetometer to a monopole. The magnetic field outside 
of the ends of a long solenoid is similar to that produced by a monopole. A solenoid can thus be 
considered as possessing two oppositely charged "pseudopoles" of pole strength g = N-I-S/go 
in units of the Dirac magnetic charge. Here iV is the number of turns per metre length, / is the 
current and S is the cross sectional area of the coil and go = 3.3 x 10~ 9 Ampere-metres is the 
Dirac magnetic charge introduced above. Hence the current and radius of the solenoid can be 
chosen to mimic the desired pole strength. 2 To calibrate, the solenoid was stepped through the 
magnetometer. Data were taken with different currents subtracting the measurements with zero 
current to correct for the dipole impurities in the coil and its former. The measured increase 

'Here z is the longitudinal coordinate defined with z = as the nominal positron-proton interaction point and 
with the positive z axis along the proton beam direction. 

2 A numerical study integrating the Biot-Savart equation for the magnetic field outside the dimensions of the 
magnetometer coil showed that this simple formula is accurate to better than ±3%. 
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in current in the magnetometer following the passage of one end of the solenoid was found to 
vary linearly with the pseudopole strength. Four coils, the details of which are given in table[H 
were used at different times for the calibration procedure. Fig. shows the magnetometer 
current divided by the solenoid pole strength as a function of the pole strength in units of go- 
There is a good consistency between the calibrations and the magnetometer is linear over more 
than 2 orders of magnitude in pole strength. The uncertainty on the point at the lowest pole 
strength is large because the current step size is small at such a low excitation of the calibration 
coil. Small differences, at the level of ~ 0. 03g£>, of the current readings occurred between 
traversals, presumably due to the system picking up specks of slightly magnetised dust between 
the traversals. Such differences show up as noise for low excitations but are less important at 
higher excitation of the calibration coil. The mean values of the calibration factor, together with 
the root mean square deviations, from each coil are given in tabled 

To simulate trapped monopole behaviour the long solenoid was placed along a beam pipe 
strip and both were passed (jointly) through the magnetometer. Only one end of the calibration 
solenoid was allowed to traverse the magnetometer hence simulating the passage of a monopole 
in the strip. Fig. |3] shows the absolute value of the measured magnetometer current as the strip 
alone was stepped through and when pseudopoles of values 2.3go and — 2.3gn were attached to 
the strip. The large structure at the centre comes from the magnetic fields from the permanent 
magnetic dipole moments in the aluminium. The final current persists when the pseudopoles are 
present. When the pseudopoles are absent the current returns to zero despite the very large per- 
manent dipole moments. The inset of Fig. |3] shows the value of the measured current (on a lin- 
ear scale) as the strip leaves the magnetometer coil. The values of current for these pseudopole 
strengths are equal and opposite and at the value expected from the calibration performed purely 
with the solenoid. This proves that a magnetic monopole attached to the beam pipe would have 
been detected by the magnetometer. 



3 Results 

3.1 Magnetometer Scans 

The data were taken in four separate sets, taking just over one day for each set: in December 
2002, January 2003, May 2003 and January 2004. In the first two sets of data (Dec. 2002 and 
Jan. 2003) all the strips from the central beam pipe, covering —0.3 < z < 0.3 m, were passed 
through the magnetometer once for each measurement. The values of the residual persistent 
current were computed from the difference between the first reading, typically ~ 20 cm before 
the strip entered the magnetometer, and the last reading which came typically ~ 30 cm after 
the strip left the magnetometer. These were then converted to Dirac Monopole units (g D ) by 
dividing by the calibration constant, determined as described above. The results for the long 
strips are shown in Fig. |4j In the first dataset (Dec. 2002) only single measurements were made 
on each long strip (except strip 13) and these are shown as open circles in Fig. |4j Two of the 
strips measured showed persistent currents of value expected from the passage of a magnetic 
charge of about +0.7gr>. Here a positive pole is defined to be a North seeking pole, i.e. one that 
is accelerated in the +z direction by the HI magnetic field. All the strips were then remeasured 
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several times in the second set of data (Jan. 2003), shown as closed circles in Fig. None 
of them (except a single reading for strip 3) showed a persistent current after traversal through 
the magnetometer. It was therefore assumed that the observed persistent currents during the 
first set had been caused by random jumps in the base line of the magnetometer electronics. 
It can be seen from Fig. |4| that none of the strips showed a persistent current which appeared 
consistently in more than one reading. Fig. |5] shows the results of the measurements on the 32 
short segments from the central beam pipe. Sample 6 showed a reading at ~ 0.5go for the first 
measurement but the measurement was compatible with zero when it was remeasured. It was 
therefore assumed that the first reading was due to a base line shift. We concluded that none 
of the long strips or short segments showed any consistent signal for a monopole in multiple 
readings and therefore that there were no magnetic monopoles trapped in the central section of 
the beam pipe. 

In the third and fourth sets of data (May 2003, Jan. 2004) the strips from the downstream 
beam pipe were investigated. In these sets of data the magnetometer measurements proved to 
be less stable. This was probably due to the induced currents from the large permanent dipole 
moments encountered which caused the magnetometer to lose its memory of the zero level. 3 
The permanent dipole moments in the downstream beam pipe were found to be much larger 
than those in the central section. These dipole moments were all observed to be aligned along 
the HI magnetic field, i.e. in the same direction as the proton beam. In these datasets the strips 
were passed through twice, first with the strip length parallel to the proton beam, termed +z 
end first, and then with the strip aligned in the opposite direction, i.e. the length antiparallel 
to the proton beam direction, termed —z end first. Fluctuations of the base level of size of 
0.7 gn were observed to happen much more frequently than for the central beam pipe section 
(Dec. 2002 and Jan. 2003 data) and in a more systematic way. Fluctuations of —0.7g D were 
consistently present for the traversal with the —z end first and fluctuations of +0.7gD were 
observed in the first 25% of the traversals with the +z end first but absent in the remaining 
75% of the readings in this orientation. To check that a monopole would have been seen, even 
in these adverse conditions, the long calibration coil was placed on a strip and the procedure 
repeated. The expected deviations from the calibration coil were seen superimposed on the 
base line shifts, confirming that a trapped monopole would have been seen had it existed, even 
in these conditions. 

The permanent dipole moments in the downstream strips gave readings in the magnetome- 
ter up to 3 orders of magnitude larger than the final persistent current expected from a Dirac 
monopole. This is larger than those seen in Fig.|3J It was found that the base line shifts could 
be avoided by demagnetising the strips in a low frequency decreasing magnetic field of initial 
strength 0.1 T. This is less than 10% of the HI magnetic field. The binding energy of monopoles 
in aluminium, the main constituent of the beam pipe, is thought to be hundreds of keV [36] com- 
pared to those of atoms which are at the eV level. Hence it was thought that such magnetic fields 
would be unable to dislodge a trapped monopole. All the strips were therefore subject to such 
a demagnetising field and remeasured. Demagnetisation was found to reduce the permanent 
dipole moments in the aluminium strips by about a factor of 20. After this procedure no further 
base line shifts were observed after passage of a strip. 

Fig.|^lshows all the readings of persistent current after demagnetisation plotted against sam- 
ple number. The strips were passed through several at a time for these data. Sample 17 consists 

3 The cause of this was not understood. 
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of the 13 strips of the central beam pipe (shown individually in Fig. |U) passed through the 
magnetometer as a bundle. It is concluded from Figs. |5] and that no monopole of strength 
greater than 0.1 Dirac magnetic charge unit had stopped in any of the measured pieces which 
constituted 93 ± 3% of the total beam pipe. The remainder was lost in the cutting procedure. 



3.2 Upper Limits on the measured Cross Sections 

To derive an upper limit on the measured cross section it is necessary to compute the acceptance, 
i.e. the fraction of the monopoles produced which would have been detected. A model of the 
production process is therefore needed. Two models were used to compute the acceptance by 
Monte Carlo technique. In each of these a monopole - antimonopole (MM) pair was assumed 
to be produced by a photon-photon interaction. The first model (model A) assumed spin 
monopole pair production by the elastic process e + p — > e + MMp through the interactions of 
a photon radiated from each of the electron and proton. The proton was assumed to have 
the simple dipole form factor 1/(1 + Q 2 /0.71GeV 2 ), where Q 2 is the negative square of the 
four momentum transferred to the proton. The second model (model B) assumed spin 1/2 
monopole pair production by the inelastic process e + p — > e + MMX (where X is any state) 
through a photon-photon fusion interaction with a photon radiated from the electron and one 
radiated from a quark in the proton. The photon is radiated with a simple distribution given 
by (1 — r)) 5 /r], with i] the fraction of the proton's energy carried by the photon. While the 
models implement the kinematic correlations in each event it should be noted that they depend 
on perturbation theory and therefore the predicted cross sections are unreliable, as mentioned 
previously. Events were generated according to model A using the programme CompHEP [38] 
and using a dedicated programme for model B. The generated final state particles were tracked 
through the HI magnetic field to the beam pipe. If the thickness of beam pipe traversed was 
greater than the calculated range of the monopole in aluminium, it was assumed to stop. In 
this way the fraction of monopoles, which were detected by stopping in the beam pipe, was 
computed. 

Monopoles experience a force gB in a magnetic field B. With the field aligned along the z 
axis they have a parabolic trajectory with 

/ x olBlr 2 r 

z (r) - z v = 0.5 y -^- + —— (1) 
ePrpr tan 6*0 

where z v is the z coordinate of the vertex and z(r) is the coordinate of a point on the trajectory 
at distance r from the proton beam. The transverse momentum and tranverse velocity of the 
monopole are Pt and /3y, respectively. The initial angle of the monopole to the proton beam 
direction is 9 and e is the unit of electric charge. In this equation g is the magnetic pole 
strength which is negative (positive) for South (North) poles which decelerate (accelerate) in 
the +z direction in the HI magnetic field. The geometric acceptance is the fraction of the 
monopoles which traverse the beam pipe in the sampled length. The total acceptance is this 
fraction times the fraction which stop in the pipe. The range of monopoles in aluminium was 
computed by integrating the stopping power, dE/dx, given in [28] adjusted for the electron 
density in aluminium. Fig.|7]shows the computed range (normalised to mass), for monopoles of 
strength g D , versus P/M = j3~f where P and M are the momentum and mass of the monopole, 
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respectively, and /3, 7 are its velocity factors. The stopping power was computed in [27-29] by 
classically considering the long range monopole interactions with atomic electrons and hence 
the result is unaffected by the limitations of the use of perturbation theory. 

Figs. [U and |9] show the total efficiency for stopping a monopole in the beampipe for models 
A and B, respectively. These were computed for magnetic charges of 1, 2, 3 and 6go, using the 
range calculations shown in Fig. 13 divided by the square of the monopole charge considered. 
The choice of magnetic charges was motivated by the Dirac quantisation condition [1] or the 
Schwinger modification [18] applied to the electron as the fundamental unit of electric charge 
(magnetic charges = lgo and 2g D ) or to the d quark as the fundamental unit of electric charge 
(magnetic charges = 3gc and 6go)- The acceptance increases rapidly as the magnetic charge 
increases since larger charges have higher dE / dx so that a greater fraction of the monopoles 
stop in the beam pipe. Hence the curve for 6go will also be approximately the acceptance 
for higher charged monopoles. The acceptance for South poles is somewhat larger than that 
for North poles since they decelerate in the HI magnetic field, losing some energy, so that 
they stop more readily in the beam pipe. Higher mass monopole pairs are produced at smaller 
angles to the proton beam and tend to hit the downstream beam pipe. However, they are more 
energetic than for lower masses. Hence high masses with low magnetic charge pass through 
the downstream pipe whereas higher magnetic charges stop. This accounts for the rise in the 
acceptance at higher masses for magnetic charges of 2g D and 3gr>. The efficiencies for model 
A tend to be smaller than those for model B since in the latter the monopoles have a smaller 
mean transverse momentum than in the former which leads to a greater fraction of monopoles 
stopping in the beam pipe. 

The upper limit on the cross section for monopole-antimonopole pair production was de- 
rived within the context of each model, as follows. The failure to observe a monopole candidate 
means that there is an upper limit of 3 monopole pair events produced at the 95% confidence 
level. The cross section upper limit is then calculated from this, taking into account the uncer- 
tainties in the measured integrated luminosity, in the fraction of the pipe surviving the cutting 
procedure, and the statistical uncertainty in the acceptance computed from the models described 
above. Here the acceptance is the fraction of the monopole pairs which produce either one or 
both monopoles which stop in the beam pipe. Fig.[TO|shows the upper limit on the cross section 
at 95% confidence level for monopoles of strength 1, 2, 3 and 6go using acceptances determined 
from model A. Fig.fTTIshows the upper limits determined using the acceptances from model B. 

Several other experiments have also produced limits on monopole production cross sections 
for different masses and charges [9, 10, 12-16, 33]. However, owing to the lack of a reliable 
field theory for monopole production, different model assumptions were made in their deriva- 
tions. Furthermore, although a universal production mechanism for monopole production can 
be postulated, comparisons of cross section limits in processes as diverse as e + p, pp and e + e~ 
are difficult. This is the first search in e + p collisions. It could provide a sensitive testing ground 
for magnetic monopoles if a monopole condensate is responsible for quark confinement [6]. 

4 Conclusions 

A direct search for magnetic monopoles produced in e + p collisions at HERA at a centre of 
mass energy of y/s = 300 GeV has been made for the first time. Monopoles trapped in the 
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beam pipe surrounding the interaction point were sought using a SQUID magnetometer which 
was sensitive down to 0.1 Dirac magnetic charges (0.1g D ). No monopole signal was observed. 
Upper limits on the monopole pair production cross section have been set for monopoles with 
magnetic charges from 1 to 6go or more and up to a mass of 140 GeV within the context of the 
models described. 
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Table 1 : Description of the calibration coils. 



Coil 


1 


2 


3 


4 


Core diam.(mm) 


6.08 


3.15 


2.1 


2.1 


Coil length (mm) 


700 


700 


300 


300 


Wire diam.(mm) (including insulation) 


0.18 


0.18 


0.1 


0.1 


Turns per metre 


11000 


11000 


10000 


30000 


Coil area (S mm 2 ) (including wire) 


32.6 


9.7 


3.80 


4.54 


Uncertainty in area 


3.3% 


6.0% 


10% 


10% 


Mean magnetometer current per go 


1.10 


1.09 


1.11 


1.11 


(arbitrary units) 










R.M.S. deviation of the readings 


0.016 


0.024 


0.030 


0.017 
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Superconducting 
Coil 




Sample 
(fixed to belt) 



Conveyor Belt 



j \_ 



SQUID 

(sensor and electronics) 



Figure 1: The schematic diagram shows the principle of the method. The conveyor belt travelled 
in steps of typically 5 cm until the sample traversed completely the superconducting coil. At 
each step the conveyor belt stopped for 1 sec before the current in the superconducting coil 
(magnetometer current) was read to avoid the effects of eddy currents. The time for each step 
was typically 3 sees. 
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Figure 2: The magnetometer current divided by the calibration pole strength as a function of 
the calibration pole strength for the four coils used (see tabledfor details of the coils). 
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Figure 3: The absolute (unsigned) value of the calibrated magnetometer current on a logarith- 
mic scale versus step position (z m ) for a strip from the central beam pipe region (—0.3 < z < 
0.3 m). The solid line shows the measurements with the long strip alone. The closed (open) 
points show the measurements with the long strip together with the calibration solenoid excited 
to simulate a pole of strength +2.3gD (—2.3go)- The inset shows the signed measurements of 
the calibrated magnetometer currents versus the step position for z m > 0.8 m on a linear scale. 
The expected persistent currents for monopoles of strength ±2.3g D are shown by the arrow on 
the logarithmic plot and by the numbers in the margin on the inset linear plot. 
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Figure 4: The measured persistent current (in units of go) in each strip, after passage through 
the magnetometer, plotted against strip number for 13 strips of the central beam pipe. Some of 
the strip numbers are offset for clarity. It can be seen that none of the fluctuations observed in 
single readings occured consistently in other readings on the same strip showing that no trapped 
monopole was present. 
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Figure 5: The measured persistent currents (in units of g D ), after passage through the magne- 
tometer, plotted against sample number for the two strips of the central beam pipe which were 
cut into short segments. It can be seen that none of the fluctuations observed in single readings 
occured consistently in other readings on the same sample showing that no trapped monopole 
was present. 
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Figure 6: The measured persistent currents in the long strips in units of go, after passage through 
the magnetometer, against sample number after the samples had been demagnetised (see text). 
Samples 1-16 consisted of several long strips (usually two or three) from the downstream beam 
pipe bundled together. Sample 17 consisted of the thirteen long strips of the central beam pipe 
bundled together. These are shown individually in Fig. [Hbefore demagnetisation. All pieces of 
beam pipe tested are included in this plot except the short segments shown in Fig. |5J None of 
the readings indicate the presence of a magnetic monopole. 
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Figure 7: The ratio of the range to mass of a monopole of charge g D in aluminium versus (3^. 
The range was calculated from the stopping power, dE/dx, in Fig. 1 of [28], adjusted to the 
electron density in aluminium. 
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Figure 8: The efficiency for stopping monopoles of strength go, 2gD, 3go and 6go or more 
computed according to the monopole pair production model A. 
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Figure 9: The efficiency for stopping monopoles of strength go, 2gD, 3go and 6go or more 
computed according to the monopole pair production model B. 
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Figure 10: Upper limits on the cross section, determined within the context of model A, for 
monopole- antimonopole pair production in e + p collisions as a function of monopole mass for 
monopoles of strength g D , 2g D , 3go and Qg D or more. 
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Figure 11: Upper limits on the cross section for monopole-antimonopole pair production in 
e + p collisions, determined within the context of model B, as a function of monopole mass for 
monopoles of strength g D , 2g D , 3go and Qg D or more. 
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